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ABSTRACT 

We present a new method to determine Thi, the total mass-to- neutral gas ratio in 
Lyman-limits systems. The method exploits the relation between the neutral hydrogen 
column density and the magnification of background sources due to the weak gravita- 
tional lensing that these systems induce. Because weak lensing does not provide a direct 
measure of mass, one must use this relation in a statistical sense to solve for the aver- 
age Thi and its distribution. We use a detailed mock catalog of quasars (sources) and 
Lyman-limit systems (lenses) to demonstrate the applicability of this approach through 
our ability to recover the parameter Thi- This mock catalog also allows us to check for 
systematics in the method and to sketch its limitations. For a universal constant Thi 
and a sample of tiqso quasars, we obtain an unbiased estimate of its value with 95% 
confidence limits (independent of its actual value) of Thi ± 140 y^C^/nQso- 

Subject headings: galaxies: formation — quasars: absorption lines — dark matter — gravitational 
lensing 



1. Introduction 

The Lyman limit systems of neutral hy- 
drogen column density Nhi > 2 x 10 17 cm -2 
are believed to be closely related to galax- 
ies at all redshifts. Among these the rare 
damped Ly—a systems (DLAS) with Nhi > 
2 x 10 20 cm -2 are usually identified with galac- 



tic disks (Wolfe et al. (1986); see, however, 
Rauch et al. (1997)), while lower column den- 
sity systems are more likely to correlate with 
galactic or proto-galactic halos (Bahcall &: 
Spitzer 1969; Churchill et al. 1996; Mo & 
Miralda-Escude 1996; Abel & Mo 1998). Ha- 
los are the basic building blocks of structure 
formation when it comes to dark matter (DM) 



gravitational collapse and virialization in any 
hierarchical scenario. Hence, low column den- 
sity Lyman limit systems (hereafter simply 
LLS) provide a unique opportunity to study 
the connection between the DM component 
and the gas component at different redshifts 
while the two are still in hydrodynamical equi- 
librium. 

The relation between the DM component 
and the gas component in LLS can be charac- 
terized in terms of the "mass-to-gas" ratio, in 
analogy to the mass-to-light ratio commonly 
used for galactic disks. In the large scale 
structure picture this ratio represents the 
mean value of the biasing factor, which con- 
nects the spatial distribution of the baryons 
(gas as well as stellar disks) with that of the 
dark matter. Generally, the biasing may be 
a linear relation between the gas over-density 
and the DM over-density, or a more elab- 
orated functional connection (e.g., Dekel & 
Lahav 1999). Biasing can be deterministic, 
but in reality it is more likely to be stochas- 
tic. It is possible to calculate such parameters 
as the mean biasing in any scheme in which 
there is a functional relation between the bias 
and the overdensity. 

In order to determine the mass-to-gas ra- 
tio, one has to observe at least two quantities, 
where one of them is directly related to the 
gas density and the other probes the DM den- 
sity field. This has been in part applied by 
Croft et al. (1998) who presented a method 
by which they recover the DM density fluc- 
tuations from the Ly— a forest. A detailed 
comparison to iV-body simulations provided 
them with the "experimental" bias parame- 
ter and spared the need to solve for it inde- 
pendently. Based on these simulations, Katz 
et al. (1996) and Hernquist et al. (1996) previ- 
ously realized that Ly— a absorption systems 
consist of a variety of physical objects: col- 
lapsing halos, cusps in phase space for density 
shells, and galactic disks. The nice agreement 



between the results from the simulations and 
the observed column density function added 
credibility to this picture. 

The appeal to iV-body simulations can be 
computationally expensive, and it is unclear 
whether the higher resolution needed for the 
LLS can be achieved today in a large cos- 
mological simulation (Gardner et al. 1999). 
The simulations also make assumptions on the 
physical conditions and the equation of state 
of the gas, which also limits their use. It is 
therefore desirable to construct an indepen- 
dent method by which the mean mass-to-gas 
ratio and other parameters can be measured. 

Here we propose to use gravitational lens- 
ing to determine the matter content, baryonic 
fraction and ionization state of the LLS. The 
absorption features, being a direct probe of 
the neutral hydrogen content in LLS, provide 
complementary data needed to reveal the con- 
nection between the matter and HI densities. 

The absorption features are well known 
and well measured (see Petitjean 1998; Rauch 
1998, for recent reviews). Their interpreta- 
tion, although subject to uncertainties in fit- 
ting the line profile, is rather straightforward 
and can be expressed in terms of their col- 
umn density Nhi, their equivalent width, and 
their velocity parameter. The lensing due to 
the moderate column densities of the LLS has 
not previously been considered because of the 
low signal expected. 

So far only lensing by DLAS (N H i > 2 x 
10 20 cm -2 ) has been studied for its effect on 
the magnification bias for the background 
source population, as well as on the distri- 
bution of the observed impact parameters 
with respect to the center of the lens (Smette 
et al. 1997; Bartelmann & Loeb 1996; Smette 
1995). Lensing by column densities of Nhi < 
10 20 cm~ 2 seemed unreasonable to consider 
because of the very small effect it has on 
the background sources; a typical quasar is 
only magnified by about one percent. The- 
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oretically, only the overall effects of lensing 
on quantities like the luminosity function of 
quasars have been considered (Hamana et al. 
2000; Pei 1995a), while observationally a pos- 
sible excess in the number counts of metal line 
systems has been detected (vanden Berk et al. 
1996). 

Other manifestations of weak lensing have 
been detected and used to map the mass dis- 
tribution in clusters and larger scale structure 
(see Bartelmann & Schneider 2000 for a re- 
view). The standard lore in weak gravita- 
tional lensing is to use the ellipticity as in- 
duced by the intervening gravitational shear 
field in order to resolve the projected mass 
distribution along the line of sight (Kaiser 
& Squires 1993). Due to the small level of 
magnification and the unavailable absolute lu- 
minosity of the sources, the observed flux of 
the background sources has never been used 
in weak lensing studies. For this reason, 
point sources, for which ellipticities cannot be 
measured, cannot be used to map interven- 
ing mass concentrations in this fashion. The 
weak magnification of point sources has usu- 
ally been considered only as a nuisance in 
extracting information from the background 
sources (e.g., in the case of supernovae la 
Wambsganss et al. 1997; Porciani & Madau 
2000). 

Until now, a key aspect of the phenomenon 
of weak lensing by LLS has been neglected, 
the fact that the same systems that cause 
the magnification are also responsible for the 
troughs in the continuum emission of the 
background source. Both are proportional to 
projected surface densities along the line of 
sight. The two are therefore correlated. It is 
this correlation, and the fact that it holds for 
each individual system that we wish to exploit 
in order to extract statistically the relation 
between the total mass of these systems and 
their HI content. In addition, the absorption 
system gives the redshift of the lens. 



In a large sample, the magnification of 
point sources can have a strong statistical ef- 
fect. Treating all quasars as a set eliminates 
in part the need to know the absolute mag- 
nitudes of the individual sources. The mag- 
nification shifts the luminosity function of the 
combined population in a calculable way. The 
lensed population can be "de-magnified" (as- 
suming a mass-to-gas scheme) on an individ- 
ual basis. The criterion for a successful correc- 
tion is the return to the statistical properties 
of an un-lensed population. 

Every statistical approach requires a large 
data base in order to maximize the signal 
to noise ratio. Fortunately in the upcoming 
years such a data base will become available 
when the quasar sample from the Sloan Dig- 
ital Sky Survey (SDSS) is completed. This 
sample should serve as an excellent target list 
by itself or for follow-up observations. Com- 
bined with high resolution spectra for a se- 
lective sub sample from ground- and space- 
based telescopes it will allow coverage of 
Ly— a absorption features over a wide range 
of redshifts. 

Here, we propose a recipe for an obser- 
vational project to determine the mass-to- 
neutral- hydrogen ratio in LLS, check different 
bias schemes, or apply a similar technique to 
the mass-to-metal-line ratio in these systems. 

In §2 we define our notation and present 
the formalism of the technique. In order to 
describe the main features of the proposed 
method in §3, we introduce a toy model which 
avoids many complications and which allows 
us to outline the ingredients of the technique. 
The next section (4) is devoted to the gener- 
ation of a realistic mock quasar and Lyman- 
limit system catalog which enables various es- 
timates of the method's biases and errors, as 
well as the functional dependence of the ob- 
tained signal-to-noise ratio on the input pa- 
rameters. In the same section, we use Monte- 
Carlo realizations of that catalog and describe 
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our procedure for measuring the mass-to-gas 
ratio. The method is implemented, and its 
advantages and limitations are discussed. We 
conclude this section by additional checks that 
verify that the measured effect indeed comes 
from gravitational lensing. In section 5 we 
sketch the preferred procedure for actual ob- 
servations, consider extensions of the method, 
and estimate additional effects like dust. Sec- 
tion 6 summarizes the method's strengths and 
weaknesses. 

2. The Theory 

It is common in gravitational lensing to ex- 
press the magnification in terms of two dimen- 
sionless quantities; re, the convergence and 7, 
the shear. The magnification [i is then given 
by 



(1 



Y 



(1) 



(Narayan &i Bartelmann 1999; Schneider et al. 
1992). In the weak lensing limit, where re ~ 
7 1 this equation becomes 



H = 1 + 2k. 



(2) 



The convergence, re, is the projected surface 
density, S, along the line of sight divided by 
the critical surface density, 



4vrG 



D c s ; Arfr = 



DiDis 



(3) 



where Di tS j s are the angular diameter dis- 
tances between observer, lens and source. 

For absorption systems the column density 
in neutral hydrogen, Nhi, along the line of 
sight is related to £ via 



Nhi=Y H i^-/J 9 £, 



(4) 



where /h the fraction of baryons that are in 
hydrogen, toh is the hydrogen mass, /& is the 
fraction of total mass made up by baryons, 



and f g is a geometry factor that relates the 
possible different geometries of the baryons 
and the dark matter (e.g., gas disk in a spher- 
ical dark matter halo, etc.). The fraction of 
neutral hydrogen is defined by 

v Phi N ffl 

y hi = = -Tj— , (5) 

PR N H 

where in the last equality we explicitly assume 
that the volume density fraction is identical 
to the surface density fraction. The geometry 
factor, f g is of order unity if LLS can be ade- 
quately described by a roughly spherical halo 
with the DM and baryons well mixed. 

We can define a projected mass-to-gas ra- 
tio, 



Thi = 



m H iN H i 



(6) 



which is determined by these different factors. 
Observations suggest values of /& ~ 0.1 (Kirk- 
man et al. 2000), fu — 0.75, and Yhi in the 
range of 10~ 2 — 10~ 3 for z = — 5 (Prochaska 
1999; Haardt & Madau 1996). Hence from 
Eqs. (4) and (6) we expect Yhi to be in the 
range of 10 2 - 10 4 . 

In a similar manner it is possible to define 
a mass-to-ionic-column-density ratio that re- 
lates the magnification to a metal line column 
density. This introduces an additional param- 
eter, the metallicity of the gas, or, more ac- 
curately, the specific abundance of a certain 
metal in a specific ionization state. Through- 
out this paper, except in §5, we will focus 
specifically on HI absorption. However, ev- 
erything we will present can also be applied 
to any observed metal line. 

With the above definitions we can now ex- 
press the magnification of a quasar by an ab- 
sorption system along its path as 

M=1 + 2 T HI N HI m H (?) 

2_/ cr 

If we could measure the magnification of a sin- 
gle quasar, we could determine Yhi for every 
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observed absorption system. However since 
there is no way to know a quasar's un-lensed 
luminosity, we must exploit the correlation be- 
tween magnification and the observable Nhi 
to determine Thi statistically. 1 

For a statistical sample of quasars, mag- 
nification has two effects. It causes the ob- 
served solid angle to correspond to a smaller 
solid angle in the source plane, and therefore 
the density of quasars with fluxes intrinsically 
above a flux limit decreases. However, magni- 
fication also increases the flux received from 
each quasar and therefore more quasars are 
detected above the flux limit. The net effect 
depends on the slope of the luminosity func- 
tion of quasars. An increase in the number 
of quasars is detected if the slope is steeper 
than —2 and a decrease otherwise. In the sit- 
uation we will be considering, magnification 
will tend to increase the number of quasars 
brought into a sample. 

In the next section, we illustrate the 
essence of the effect studied here with a toy 
model. Most of its simplifications will be 
dropped in later sections, to which readers fa- 
miliar with gravitational lensing might want 
to skip right away. 

3. Simple Model 

In this section we introduce a simple model 
that demonstrates the basic concepts of the 
method for a statistical recovery of Thi with- 
out the complications that will be introduced 
in the following, more realistic case. 

Let us assume all quasars are at z = 3 and 
their luminosity function (LF) is given by a 
power law of slope —3.0. The number of ab- 
sorption systems as a function of column den- 
sity and redshift with column densities in the 

1 Strictly speaking we determine the cross-section 
weighted average of Thi. 



range 18 < logN H i < 19. 7 2 is taken to be 

dra = 0.096N^ 8 L5 (1 + z) L55 dNi 8 dz (8) 

where N±s = NHi/(10 18 cm~ 2 ). The choice of 
these parameters will be justified in section 
4.1. Throughout this paper we adopt the cos- 
mological parameters Ho = 65kms _1 Mpc _1 , 
Q.M = 0.3 and Sl\ = 0.7. 

Using these distribution functions, we con- 
struct a sample of quasars and the associ- 
ated absorption line systems along their lines 
of sight. The magnification value for each 
quasar is calculated according to the systems 
in its foreground, and the apparent luminos- 
ity of the quasar is correspondingly corrected 
by applying Eq. 7. The dilution of quasars 
by solid-angle magnification is taken into ac- 
count by randomly removing quasars. This is 
performed for values of Thi of 0, 1000, and 
2000 (we loosely denote the un-lensed case by 
Thi = 0). 

A flux cutoff at 45 /uJy (B brighter than 
20 m ) is then applied to the magnified sample. 
Quasars which initially were below this flux 
limit are magnified by their LLS, exceed the 
flux limit, and introduce more high column 
density systems into the sample. 

Figure 1 shows the quasar LF in this toy 
model. The upper lines depict the LF of all 
quasars in the sample with the three Thi val- 
ues. Since the fraction of significantly mag- 
nified quasars is negligible, the magnification 
(different Thi values) has no detectable effect 
on the total LF, in accord with the conclu- 
sions of Pei (1995a). The lower three lines 
depict the LF as derived from a subset of 
quasars with at least one absorption system 
with log Nhi > 19.0. For this subset it is ev- 
ident that the LF is shifted to the right with 
increasing Thi, he. increasing magnification. 

2 Throughout this paper logarithmic column densities 
are assumed to be measured in cm~ 2 . 
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Fig. 1. — The observed luminosity function 
of our toy model for all quasars and for 
quasars with at least one absorption system 
with logNni > 19.0. Each is plotted for three 
values of T m (0, 1000, 2000). The upper three 
lines lie on top of each other as weak lensing 
has a negligible effect on the total quasar lu- 
minosity function, as was previously found. 
However the magnification by gravitational 
lensing is discernible in the subset of quasars 
that contain high column density absorption 
systems. 

In principle, the relative heights of the dif- 
ferent lines in Fig. 1 are proportional to Thi 
and could be used to find its value. Unfortu- 
nately, in practice there is only one realization 
of the Universe so only one of the lines (for 
each Nhi cutoff) in Fig. 1 is observed. In or- 
der to find out what is the relative amplitude, 
the lowest line (un-lensed, Thi = 0) should be 
inferred from the expected fraction of high col- 
umn density systems, based on lower column 
density systems and Eq. (8). Thus the ef- 
fect of magnification is to change the fraction 
of observed quasars with high column density 
absorption systems. 

This effect is more clearly seen if instead 
of looking at the luminosity function we look 
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Fig. 2. — The distribution of column densities 
of observed absorption systems for values of 
T H i of 0, 1000, and 2000. The number of ab- 
sorption systems observed per 1000 quasars is 
plotted versus the log of the column density. 
For higher values of Thi more high column 
density systems are brought into the sample 
because of the magnification bias on the back- 
ground sources. 

at the distribution in column density of the 
absorption systems. Figure 2 shows this dis- 
tribution for different values of Thi- The ef- 
fect of increasing Thi is immediately obvi- 
ous; more high column density systems are 
brought above the flux limit due to magnifica- 
tion bias. The distribution also deviates from 
the original power law as the effect becomes 
progressively more pronounced with increas- 
ing N H i. 

The magnification which causes the line di- 
vergence in Fig. 2 is not only a function of 
Thi but also of the lens' redshift through D e ^ 
in Eq. (3). We demonstrate this dependence 
in Figure 3 where the probability of encoun- 
tering an absorption system (here normalized 
to integrate to unity) as given by Eq. (8) is 
weighted by the redshift dependent factors of 
Eq. (3) that contribute to the magnification. 
The resultant function is a modification of Eq. 



6 



0.8 
£> 0.6 



1 1 1 1 1 1 1 1 
_ / " \ 


| 1 1 1 1 | 1 1 1 1 | 1 1 1 1 | 1 1 1 1 

weighting function 

weightedP abs 


/ 

/ 

/ 

/ / 


\ 

\ \ 
\ \ 
\ \ 


17 

\ 


X \ 

s \ 

x \ — 

X \ 

x \ - 

X \ 



0.0 0.5 1.0 1.5 2.0 2.5 3.0 



z 

Fig. 3. — Magnification weighted probability 
of an absorption system ciS cl function of ab- 
sorption redshift (thick solid line) for quasars 
at z = 3. The unweighted probability of an 
absorption system (Eq. 8, thin dotted line is 
weighted by D c r of Eq. (3) (dashed line). All 
three functions are normalized so that they 
integrate to unity. The combination of mag- 
nification weighting which peaks at interme- 
diate redshifts, and a monotonically increas- 
ing probability of encountering an absorption 
system with redshift, results in a rather broad 
redshift range in which lensing by absorption 
systems is important. 

(8) (again normalized to integrate to unity) 
with a broad distribution in redshift. This 
means that lensing by absorption systems is 
relevant over a wide range of redshifts. 

This redshift dependence distinctly affects 
the counts of absorption systems per unit red- 
shift for a fixed number of quasars. This effect 
depends on the value of the magnification co- 
efficient, Thi- Figure 4 shows the number of 
absorption systems with logNni > 19.0 as a 
function of absorption redshift. The increase 
in number of high column density absorption 
systems happens preferentially between red- 
shifts of 0.5 and 2.0. This is a clear sign that 
the increase is due to lensing for which the 
magnification becomes very small as the lens 




1+z 



Fig. 4. — The distribution of observed absorp- 
tions systems with column densities greater 
then 10 19 cm -2 in redshift for Thi values of 0, 
1000 and 2000 for quasars at z = 3. Since 
lensing depends on the redshift of the lens, 
systems that are at redshifts between 0.5 and 
2.0 are preferentially brought into the sample. 

approaches the quasar redshift. A deviation 
from the redshift power-law dependence of Eq. 
(8) is thus a statistical tool by which the am- 
plitude of the magnification can be inferred. 
These effects do depend to some extent on cos- 
mology. 

The toy model we have used here highlights 
the main effects of treating the magnification 
by intervening absorption line systems in a 
statistical manner. The luminosity function 
of quasars is shifted to higher luminosities, 
a shift that can be observed if a cut by col- 
umn density along their line of sight is per- 
formed. However, it is cumbersome to get 
a reference line for a non-lensed luminosity 
function, whereas it is more straightforward 
to detect the change in the column density 
function and its redshift distribution. 
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4. Recovery method for Thi 

The results of the toy model as presented 
in the last section portrayed the basic effects 
of lensing on the statistical distributions of 
quasars and their absorption systems. In real- 
ity, however, there are various deviations from 
the idyllic picture of the toy model. In order 
to assess our ability to use these altered dis- 
tributions for the recovery of Thi we need to 
take all of these complications into account. 
Only a detailed imitation of a realistic sit- 
uation will allow us to draw conclusions re- 
garding the success or failure of the recovery 
method. 

The main complications that should be in- 
cluded in a more realistic picture are; 

• The observed luminosity function of 
quasars resembles a double power law, 
rather than a single one, and thus an ad- 
ditional scale (the "knee") is added to 
the problem (Boyle et al. 1988; Hartwick 
& Schade 1990; Warren et al. 1994; 
Boyle et al. 2001). 

• The observed quasars span a large range 
of redshifts. This together with the 
redshift dependence of the LF "knee" 
results in a complicated observed flux 
function. 

• An observed sample is flux limited, but 
quasars may be magnified above the flux 
limit depending on their absorption sys- 
tems and the value of Thi- 

• Frequently, more than one absorption 
system intersects the line of sight to a 
quasar. These events, especially for high 
column densities, contribute a substan- 
tial signal to the shifts of the distri- 
bution functions. Multiple magnifica- 
tion events must thus be included in the 
analysis. 



• Observational errors in Nhi and flux 
measurements introduce various biases 
to the method. These biases must be 
treated to evaluate their systematic ef- 
fects on the recovery result. 

For all these reasons one must appeal to 
Monte-Carlo realizations of a realistic mock 
catalog of quasars and absorption systems. 
The recovery method is demonstrated and 
tested using these Monte-Carlo realizations. 

4.1. The Mock Catalog 

We generate a sample of quasars according 
to their observed flux and redshift distribution 
dnQso/d-Sdz. Although this is the quantity 
which is directly observed, it is not reported 
in the literature. Instead it is converted to a 
luminosity function by binning under an as- 
sumed cosmology. We thus start with a fit to 
the luminosity function and recover the ob- 
served flux function. We start with the dou- 
ble power-law fit to the luminosity function as 
given by Pei (1995b). 



dn QS Q _ _ _ 

dLdV ~ (L/L z )& + (L/L z )fo 



(9) 



where L z = L* exp [— (z — z if ) 2 /{2a 2 )} de- 
scribes the redshift evolution of the lumi- 
nosity function, = 234 (comoving)Gpc~ 3 , 
Pi = 1.83, fo = 3.70, = 2.63 x 1O 13 L , 
= 2.77 and a = 0.91. We convert this to 
the observed flux function by 



dn Q so 
dSdz 



47T- 



dn 



QSO 



dL dV 



dLdV dS dz 



(10) 



For the probability of encountering an ab- 
sorption system P(Nhi, z) we assume this 
function is separable in Nhi and z as observa- 
tions seem to imply. Storrie-Lombardi et al. 
(1994) found that the number density of LLS 
increases as (l + z) L55 . Many authors (Tytler 
1987; Press & Rybicki 1993; Hu et al. 1995; 
Lu et al. 1996; Kirkman & Tytler 1997; Kim 
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et al. 1997) showed that the distribution in 
column density is well fit by Nhi -1 ' 5 over ten 
decades. Thus we take the probability density 
to be 

dP = 0.096Nf 8 L5 (l + z) L55 dNi 8 dz (11) 

with Nis = NHi/(10 18 cm~ 2 ) as in section 
3. We only consider absorption systems with 
column densities between 10 18 cm -2 and 5 x 
10 19 cm~ 2 . At the lower limit the magnifica- 
tion in the most favorable configuration (max- 
imal Z?efr and high Thi) is only 0.01 so we can 
safely neglect lower column density systems 
when considering lensing. 

The upper limit is set because we would 
like all the systems we consider to contain gas 
in a similar ionization state (Yhi), i.e, similar 
Thi- This upper limit must be set lower than 
the column density of damped Ly— a (Nhi > 
2 x 10 20 cm~ 2 ) that consist of mostly neutral 
gas (Prochaska & Wolfe 1996). A kinematic 
analysis of damped Ly— a by Mailer et al. 
(2001) suggests that cold gas disks at z ~ 3 
extend down to Nhi ~ 4x 10 19 cm -2 , and local 
observations of the average surface densities 
of gas disks yield Nhi ^ {few) x 10 19 (Zwaan 
et al. 1997). Taking inclination effects into 
account this implies that absorption systems 
with column densities below 5 x 10 19 cm -2 are 
dominated by hot gas. We therefore take 
this value as a conservative upper limit. This 
choice of an upper limit also safely keeps us 
away from strong lensing events that would 
not be described correctly by Eq. (7). 

Our results for Thi turn out to be quite 
insensitive to a change in the assumed cos- 
mology, the quasar flux function, or the less 
constrained column density function. Their 
primary effect is a slight modification in the 
number of quasars required to reach a result 
with the same level of confidence. In light 
of the current uncertainty of over an order of 
magnitude in the value of Thi, we will ignore 
these possible modifications. 



Armed with these relations we turn to gen- 
erate a mock catalog of quasars and their 
associated absorption systems. For the rare 
high (though still weak) magnification events 
which lift quasars above the flux limit of the 
survey, there is a need for a large sample 
of quasars. We have generated a sample of 
8 x 10 6 quasars and assigned them fluxes and 
redshifts according to the flux/redshift distri- 
bution, dnqso/dSdz (Eqs. 9, 10). 

In generating this catalog we have in mind 
the specifications of the SDSS and apply a 
flux limit of 45/uJy to approximate the SDSS 
limit of 19 m .7 in the g' band. The target for 
the photometric errors of the Sloan survey are 
2%. Based on what is commonly quoted in the 
literature we estimate that the errors in mea- 
suring Nhi to be around 15%. This catalog is 
used in the next section to exercise and test 
our proposed method for the Thi recovery. 

4.2. The Recovery Procedure 

Our goal is to present an unbiased recov- 
ery method for Thi, assumed at first to be 
a universal constant. Relaxation of this as- 
sumption will be discussed in §4.3. The re- 
covery is judged by the confidence levels that 
can be put on the value of Thi as function of 
its value and the number of quasars above a 
flux limit that are needed in order to get these 
levels. A magnification signal is manifested in 
the change of the distribution functions of the 
column density: the relative number density 
of different column density systems and their 
redshift distribution. For P(Nhi, z) there are 
three independent variables: the amplitude of 
the function, the slope of the separable Nhi 
part of the function, and the slope of the sepa- 
rable z part of it. Note that the amplitude can 
only be used once. We assume that the true 
underlying slopes can be read off the low col- 
umn density systems that produce only neg- 
ligible magnification, and then extrapolated 
to higher column densities. Justifications for 
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l°gio N ffl 

Fig. 5. — The distribution of the total pro- 
jected HI column density, N^ : , for values of 
T H i of 0, 1000 and 2000. Without lensing 
the distribution would be that of the Thi = 
case. Because of the magnification bias, sys- 
tems with high values of Njjj are magnified 
over the flux limit. The error bars are the 
Poisson errors from a Monte-Carlo realization 
with riQso = 10 5 . 

this assumption and ways to check its validity 
will be presented in §4.4. 

4-2.1. The total column density along the 
line of sight 

Let us now introduce a new variable, N^j, 
the total column density of HI along the line 
of sight. For the j'-th quasar in the sample, 



(12) 



where the sum is over all systems along the 
line of sight. Note that this is not weighted 
by D c g. The probability distribution of this 
variable, -P(Njjj), folds in the redshift /flux 
distribution of the quasars, the column den- 
sity /redshift distribution of the Lyman- limit 
systems, and in addition the Poisson process 
of multiple systems along the line of sight. 

The total magnification experienced by a 
quasar is the product of the magnification by 



each absorption system along its line of sight. 
Thus Ngj has a stronger correlation with mag- 
nification than the column densities of indi- 
vidual absorption systems. Another advan- 
tage of the Npjj quantity is that its measure- 
ment does not require assumptions on cosmol- 
ogy. Figure 5 shows P(N^j) for three values of 
Thi- There is a dramatic dependence on Thi 
in the distribution -P(Njjj) for higher total col- 
umn densities. This should be compared to 
Fig. 2 where the effect of lensing was much 
less evident. 

4-2.2. Rejection of the null hypothesis 

We first examine the rejection levels of the 
null hypothesis that the sample under consid- 
eration is not magnified at all. Under this 
assumption any deviation from the expected 
P(Nhi,-z) is purely because of Poisson errors 
due to a finite sample and observational er- 
rors. As discussed above the relevant dis- 
tribution function to maximize the lensing 
detection is P(N^j) rather than P(Nhi,-z). 
The equivalent projection along the z axis, 
-Pn H i( z ) = / dNHi-P(NHi, z) can in principle 
also be used, as we demonstrated in Fig. 2. 
Let us first phrase the statistics we use in the 
conventional \ 2 terms and then modify it to 
get a statistic which is more appropriate for 
the problem under consideration. 

For a sample of 71qso quasars, let n™ odel (N^j) 
be the number of expected absorber in bin j 
spanning the range [N^Nj^j + ANj^J. This 
number depends, of course, on the redshift 
distribution of the quasars in the sample. 
Let n™ odel (2;) be the number of expected ab- 
sorbers in bin i that spans the redshift range 
[z, z + Az]. If the two observed values of the 
same quantities are n° bs (N^j) and n° bs (z) we 
can construct the two x 2 functions 
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error, j 



where the sums are over the relevant bins. 
The two error contributions are from the Pois- 
son error due to the finite sample, e 



poisson 



and 

the propagating error due to the observational 
uncertainty in A^hi- The slope in the P(iVHi) 
function is affected by a Malmquist-like bias 
since it is more likely for lower column density 
systems to "leak" into a bin of higher column 
density than the other way around. We there- 
fore adopted the following form for the error 
terms, 



c poisson,j — y n j J 
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where 5Nhi is the measurement error in Nhi 
(15%). We approximate the Poisson error 
from the modeled (expected) rij since ulti- 
mately we would like to get a probability dis- 
tribution about the true underlying value of 
Thi. 

The two functions, x\ 2 are n °t indepen- 
dent as written in Eq. (13). If, for instance 
the total number of observed absorption sys- 
tems exceeds the expected one, this will affect 
both the x 2 values. In order to regain the sta- 
tistical independence, one of the x 2 functions 
should be renormalized such that the modified 
function is read 
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where the original errors, e, are re-scaled ac- 
cordingly to e. The final x 2 function is then 
expressed as the weighted sum of the two in- 
dependent functions 

x 2 = C1X1 + C 2X2 ; ci + c 2 = 1 . (16) 

In practice we find that the statistics is dom- 
inated by xi an d little improvement in the 



statistical rejection levels or confidence levels 
is obtained by including x\- We therefore set 
ci = 1 and C2 = for the rest of our analysis. 
For example, the x 2 statistics allows us to re- 
ject the null hypothesis with riQso = 2 x 10 4 
and values of Thi of 300, 400 and 500 at the 
70%, 82% and 90% confidence levels respec- 
tively. 

So far we have used standard x 2 statistics. 
However, in the current case, where we try to 
detect a non-random process (or to reject the 
absence thereof) there is a better, more ro- 
bust way of doing it. Based on Fig. 5 we no- 
tice that the high column density end of the 
Njjj distribution can be well fit by a power 
law. This power-law description is excellent 
for the non-magnified case, but it provides a 
reasonable description for the magnified cases 
as well. A full description of the fit to the 
non-lensed case is provided by its normaliza- 
tion, ao, and its slope a\. Recall that since 
is the (plain, unweighted) sum of the Nhi 
of all absorption systems on the line of sight 
to a quasar, a\ is not identical to the slope of 
dP/dN ffl of Eq. (8). The function P(Ng r ) is 
established using 10 6 un-lensed quasars from 
the mock catalog. From this we determine 
n™ odel (Njijj) for a given value of nQso and bin 
size. The best fit to the values of n™ dcl (Ngj) 
and the errors on these values as given by 
equation (14) provides the value of (ao,ai). 

In the presence of magnification, we expect 
a coherent deviation from the non-magnified 
line (ao,a>i) of the N^j distribution. This al- 
tered distribution can be fit with the two pa- 
rameters which describe it best, (oQ,a'j), by 
a standard fitting procedure, which also takes 
into account an error covariance matrix Aa 
which is a function of n™ odel (NHi) and of the 
e jS . 

At this point the new statistic we estab- 
lish compares two lines; the model line (ao, ai) 
and the fit to the observed line (a' ,a[). We 
use Monte-Carlo runs to map the probability 
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Fig. 6. — Rejection levels for the null hypoth- 
esis that the observed distribution of is 
unbiased by gravitational lensing func- 
tion of the number of quasars in the sample 
and the value of Thi- raQSO is the total num- 
ber of quasar in the sample, but only a subset 
of these quasars are used in order to calculate 
the rejection levels. The subset is determined 
by a cut in Ngj so that 19.0 < logNgj < 19.7. 
This amounts to 10% of the uqso in the case 
of no lensing and increases to 12% if Thi is 
1000. 

space in the parameters (aQ,a\) of the non- 
magnified Njjj distribution and to evaluate 
the probability of fitting a line (in log) of nor- 
malization a' Q and slope a[ given a random, 
discrete realization of this distribution with 
raQso lines of sight. The outcome is a full 
map of the distribution dP / (daodai) which is 
similar to Eqs. (13) replacing n™ odel (NHi) by 
the fit to the modeled and n° bs (N m ) by 
the fit to the realization of the NSt distribu- 
tion. Under this change of variables there is 
no longer summation over bins and the func- 
tion which effectively replaces x 2 becomes, 

= |(ao,Qi) ~ (QqXi)! 2 (17) 



Fig. 7. — Same as Fig. 6, but rejection levels 
are drawn based only on quasars with 19.5 < 
logNyj < 19.7 which corresponds to 2.1% and 
2.7% of riQso for T ffl = and T m = 1000 
respectively. 

For the z distribution this statistic cannot 
be used because magnification causes an in- 
crease of the number of systems in some bins 
and a decrease in others. Thus a magnified 
distribution is no longer well fit by a line (in 
the log). To compare instead an observed z 
distribution with a non-magnified one, we use 
the standard Kolmogorov-Smirnov test. 

Based on this ifi statistic we can draw re- 
jection levels of the null hypothesis. These are 
more stringent than those obtained using \ 2 \ 
for example with a sample of 2 x 10 quasars 
an Thi of 600 can be rejected with 99.4% con- 
fidence using the ip statistic, but with only 
96% confidence using \ 2 ■ Figure 6 shows the 
rejection levels as function of nQso and the 
true modeled value of the constant Thi- Even 
though tt-qso is of the order 10 4 — 10 5 [x axis), 
the rejection levels are drawn from a subset 
of quasars that were selected based on their 
Njijj value. This subset amounts to 10% of the 
quasars if 19 < logN^j < 19.7 and Thi = 
and increases to 12% for Thi = 1000. As ex- 
pected, the rejection confidence increases with 
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increasing nQso and Thi- For example, it 
is enough to have a total sample of 2 x 10 4 
quasars in order to reject the null hypothesis 
by more than 99% for a true Thi value of 600. 
Recall that only about one tenth of this sam- 
ple is actually being used in this calculation. 

The fraction of the sample that is actually 
being used can be further reduced by a fac- 
tor of 4 — 5 if a higher lower limit for N^j is 
adopted. Figure 7 shows the rejection levels 
if only quasars with 19.5 < logN^j < 19.7 
are taken into account. This reduces the sam- 
ple to ~ 2.1% and ~ 2.7% of the size of the 
original sample for Thi = and Thi = 1000 
respectively. However it now takes about 
4.5 x 10 4 quasars in the original sample in or- 
der to reject the null hypothesis at the 99% 
confidence level for a true Thi value of 600. 

4-2.3. Confidence levels for recovered Thi 

Rejection of the null hypothesis is only the 
first step in establishing the value of Thi and 
its functional shape. Similarly to the choice 
of C2 = in the \ 2 statistics of the previous 
section (§4.2.2, Eq. 16), we only consider the 
change in the number counts in N^j bins since 
the information in the change of the redshift 
distribution does not improve the confidence 
levels by more than a few percent 3 . 

To get the best value of Thi we de-magnify 
each quasar in the sample individually, based 
on the actual absorption systems along its line 
of sight. The observed flux, Sj of a quasar j 
at Zj is shifted to the de-magnified flux, Sj by 

/ ^ j 

UiHA z j^ z i^i ; t hi) 

where fiji is the magnification due to the ab- 
sorption system i of column density Nm* at Zj 

3 In addition to the change in the absorbers' redshift 
distribution, there is also a small shift in the redshift 
distribution of magnified (high N^j •) quasars relative 
to the non-magnified quasars. 



on the line of sight to the quasar, according to 
Eq. (7). If Sj turns out to be lower than the 
"survey" flux limit, then quasar j along with 
all its associated absorption systems are re- 
moved from the sample. The remaining sam- 
ple contains ?i<qso(^ n QSo) quasars with a set 
of n'qgQ values of N^ T . Then in order to ac- 
count for the dilution of quasars due to solid- 
angle magnification, we add quasars in the 
following way. Quasars are drawn randomly 
from the remaining sample with a probability 
1/fi and then added to the sample such that 
there are duplicates of some of the quasars 
and their N^j values. These values are fit 
(in the log) to a linear function character- 
ized by (a^a^) and compared to the non- 
magnified (i.e., the extrapolation from the low 
Njjj range) distribution fit. This is done in the 
same manner described in §4.2.2. The best 
Thi value is obtained by minimization of the 
V> statistic. 

Figure 8 shows the success of the recovery 
of Thi values along with its confidence lev- 
els. The ratio of the recovered value (T HI ) 
and the true pre-assigned value (T HI ) is shown 
as a function of T HI . Only quasars with 
19.0 < logN^Jj < 19.7 were used to draw the 
confidence levels. As the total h-qso increases 
(from the top panel down) the confidence lev- 
els of the recovery for a given true T HI value 
become narrower about this value. Note that 
the average recovered value is not biased even 
when the confidence levels are broad. 

All quantities are calculated using between 
10 — 20 Monte-Carlo runs depending on the 
value of jt-qso- With a sample of 10 5 quasars 
it is possible to measure Thi to a higher ac- 
curacy than 30% of its value if it is greater 
than 500. This is a remarkable improvement 
over the current uncertainty of an order of 
magnitude. The uncertainty in the recovered 
value of the mass-to-gas ratio does not depend 
on T^jj but only on riQso (the relative uncer- 
tainty as shown in Fig. 8 goes down because 
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Fig. 8. — The ratio of the recovered (T HI ) 
to the true (Tjjj) value of the mass-to- gas ra- 
tio and the 95% (shaded region) and 99.7% 
(dashed line) confidence regions are shown 
versus T HI for three quasar sample sizes. 



T HI increases). The 95% uncertainty in the 
recovered value of Thi is approximately equal 
to 140 A /^E. 

V "QSO 

Using only quasars with 19.5 < logNjjj < 
19.7 out of a sample size of 10 5 quasars has 
the effect of broadening the confidence levels 
to the degree of confidence level widths close 
to the case of nQso = 6 x 10 4 with 19.0 < 
log Ngj < 19.7. 

4.3. Non-constant Thi 

The assumption of a universal mass-to-gas 
ratio is quite simplistic. We here explore 
two deviations from this assumption. The 
first scheme we investigate assumes that Thi 
has a Gaussian distribution with mean Thi 
and standard deviation, cry m . The middle 
panel of Figure 9 shows the recovered value 
for Thi and the confidence levels about it for 
cjt hi = 0.5T HI . Introduction of this signif- 
icant scatter in the mass-to-gas ratio has a 



Fig. 9. — The ratio of the recovered (T HI ) to 
the true (T HI ) value of the mass-to-gas ratio 
for different mass-to- gas ratio schemes. The 
top panel is the universal constant scheme ex- 
plored in §4.2. The middle panel shows the 
result from a random mass-to-gas ratio with 
a mean of T HI and a ct hi = ^hp The bot- 
tom panel shows the results from the redshift 
dependent scheme explained in the text. For 
all three panels uqso = 10 5 and quasars with 
19.0 < logN^j < 19.7 are used with the ip 
statistic. The shaded region and dashed line 
denote 95% and 99.7% confidence levels re- 
spectively. 

modest effect on the recovery of Thi- The 
recovered mean value is slightly higher than 
T^, but this small bias is less than 50% of the 
uncertainty of the recovered value of Thi even 
for values of o"x HI as large as T HI . Thus it is 
only noticeable for small values of T HI as can 
be seen in Fig. 9. Part of this bias arises be- 
cause Thi must be positive, which skews the 
Gaussian distribution, thereby increasing the 
mean value of Thi- However, there is also the 
effect that in an observed sample higher values 
of Thi w iH preferentially be brought into the 
sample, increasing the mean Thi value of ob- 
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served absorption systems. The scatter in the 
recovered value shows no noticeable increase 
for o"x HI = 0.5T HI because the Poisson scat- 
ter in the distribution of absorption systems 
dominates the noise. Thus the statistical na- 
ture of the recovery and the large data set be- 
ing used makes random scatter unimportant 
in the recovery of Thi- 

The second scheme for the deviation of Thi 
from a universal constant is motivated by the 
ionization evolution of the Universe. If the 
Yhi parameter of LLS in Eq. (4) is a function 
of redshift then so is Thi- For the functional 
form of Yhi we adopt (Haardt & Madau 1996) 



Yhi 



Thi - Pna>m(T) 



(19) 



where Thi, the photo-ionization rate per hy- 
drogen atom, is given by 



A(l + z) B exp 



- z c ) 



21 



(20) 



with the fitting parameters A = 7 x 10 13 , B = 
0.8, s = 1.95, and z c = 2.2. am is assumed to 
be 2.6 x 10~ 13 , its value for a 10 4 K, optically 
thin gas, with density of /?h = 4 x 10~ 3 cm~ 3 . 
This roughly corresponds to the virial density 
and fb = 0.1. In practice we adjust A so that 
the maximal Thi value we obtain is 3000. 

First, we address the issue of recovering 
Thi under the assumption that it is a uni- 
versal constant. In this case we will refer to 
the true value of the mass-to-gas ratio T HI as 
the average value of Thi- For the parame- 
ters given above the average value is 42% the 
maximal value. 

When Thi really does vary with redshift, 
the bottom panel of Fig. 9 shows the results 
of the recovery method. The recovered value 
is about 2/3 the average value which reflects 
the sensitivity of lensing to redshift as shown 
in Fig. 4. Also the uncertainty in the re- 
covered value is much wider then in either the 
universal constant or random scatter schemes. 
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Fig. 10. — The number of absorption systems 
with logN m > 19.0 function of redshift. 
The thick solid line is the un-lensed case. The 
dotted line is the lensed sample magnified by 
a mass-to-gas ratio that depends on redshift 
(see text). The dot-dot-dot-dash line is ob- 
tained when the sample is erroneously de- 
magnified by an assumed universal constant 
value for Thi- De-magnifying by the correct 
z dependent scheme results in the dashed line. 

The second issue of concern is whether we 
can determine z c in Eq. (20) using our recov- 
ery method. First we need to know that Thi 
is a function of z and not a universal con- 
stant. This can be done by examining the 
distribution in redshift of absorption systems 
with logNHi > 19.0 in the demagnified sam- 
ple. Figure 10 shows that while the demag- 
nified sample fits the unbiased distribution of 
Nhi the z distribution of its absorption sys- 
tems is markedly different. 

Finding that the assumption of a univer- 
sal constant Thi provides a poor fit to the 
expected z distribution of absorption systems 
with log Nhi > 19.0, we can try other mass- 
to-gas schemes. If we assume a function of 
the form Eq. (20) we can allow z c and the 
maximum value of Thi to vary and compare 
the resulting de-magnified samples against the 
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Fig. 11. — The ability to recover the values 
of z c and the maximum value of Thi using 
the NSt distribution and the z distribution of 
absorption systems. The broad, closed region 
shows the constraints from considering the z 
distribution and the narrow arc comes from 
considerations of the distribution. This 
is done with tiqso = 10 5 . Both tests individ- 
ually allow for a large range in z c but their 
combination can yield a fairly good measure 
of these parameters. 

expected Njjj and z distributions. The results 
are shown in Figure 11 where the rp statistic 
has been used to compare to the N^j while a 
K-S test is used to determine the probability 
that the demagnified z distribution is compat- 
ible with the expected one. The narrow curve 
is the constraint from the ip statistic while the 
broader, closed region shows the constraints 
from matching the z distribution. The black 
diamond shows the true values used in the 
magnification. Values of z c that are compat- 
ible with both tests at better then the 90% 
level lie between redshifts of 1.5 and 2.9. Thus 
constraints on the redshift dependence of Thi 
can be placed by considering the z distribu- 
tion of absorbers in combination with the Njjj 
distribution. 



4.4. Additional checks for the method 

Checks that the signal measured and as- 
cribed to gravitational lensing is really due 
to gravitational lensing come naturally in our 
scheme. The first verification is to check that 
with the recovered value of Thi the demag- 
nified sample is in all ways unbiased. The 
second is to break the observed sample into 
a part favorable to being lensed and an unfa- 
vorable part and then perform the above pro- 
cedure comparing these two parts. 

The correct value of Thi should not only 
bring the distribution of N^j into agreement 
with what is expected for an un-lensed pop- 
ulation, but also the flux function of these 
demagnified quasars and the redshift distri- 
bution of their absorption systems should be 
in agreement with those quasars that are not 
significantly lensed. The ability of one value 
of Thi to bring the quasar sample into agree- 
ment with these three independent distribu- 
tions seems quite unlikely if the measured dif- 
ference in Njtt was due to some other effect 
than lensing. If the z distribution does not 
match the un-lensed one this does not rule out 
lensing, it simply rules out the dependence of 
Thi on redshift that we have assumed. A dif- 
ferent dependence that is compatible with the 
three distributions must then be found as has 
been demonstrated above. 

A second check is to divide the sample into 
those absorber-quasar pairs that are favor- 
able and unfavorable to lensing. Equation 
3 shows that lenses of the same mass with 
redshifts near the redshift of the quasar pro- 
duce less magnification. Thus separating the 
sample by this criterion creates a subsample 
that is in every way identical to the total sam- 
ple except for this one relationship that should 
only be relevant to gravitational lensing. The 
NSt distributions of the two samples can be 
compared and Thi recovered in exactly the 
same method as explained above. The advan- 
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tage of this test is that it can be made solely 
with the data at hand with no need to ex- 
trapolate P(Nhi,z) from lower column den- 
sities; however, because one divides the ob- 
served sample into two parts the signal will 
be reduced. Thus our preference is to assume 
that we can extrapolate P(Nm,z), and then 
perform the separation of the sample as an ad- 
ditional check that our assumption is correct. 
In practice one would first want to validate 
the assumption that P(Nhi, z) can be extrap- 
olated to higher column densities by check- 
ing that those systems unfavorable to lensing 
are consistent with the lower column density 
P(N ffl ,z). 

5. Observational procedure, metal lines, 
and dust 

The basis of any future large quasar study 
will be based on the five band photometry of 
10 6 quasars and spectra of 10 5 quasars ob- 
tained by the SDSS 4 . The Sloan spectra will 
identify Ly— a , Mgll and CIV lines in the 
redshift ranges, 2.3 < z < 6.4, 0.5 < z < 2.2, 
and 1.6 < z < 4.8, respectively. 

The detection of Ly— a from the ground is 
restricted to z > 2 so only about 15% of the 
absorption systems we have been considering 
in §4 will be detected in Ly— a in the Sloan 
data. Furthermore these systems will be less 
favorable to lensing because of the proxim- 
ity in redshift of the quasar and the absorber. 
However, as shown in Fig. 6, if the true mass- 
to-gas ratio exceeds 800 then with only 5000 
quasars it is possible to reject Thi = (no 
lensing) at the 90% confidence level. So it is 
possible that a lensing signal may be detected 
with just the LLS identified in the Sloan sur- 
vey. 

Mgll, CIV and other metal lines are de- 
tectable over a wide range of redshifts from 
the ground and thus from the Sloan data alone 

4 http://www. astro. princeton.edu/PBOOK/welcome. htm 



values for T]yigii and Tciv can be determined. 
While the theoretical interpretation of these 
values is more uncertain than Thi because of 
the added complication of metallicity of the 
gas, nonetheless metal lines can establish that 
there is lensing by absorption systems. The 
mass-to-metal-line ratios of different elements 
can be used to study the enrichment history 
of LLS. Using correlations between metal line 
and Ly— a column density, Thi can also be 
inferred, though with a bigger scatter. 

The optimal sample consists of LLS be- 
low redshift 2. The Cosmic Origins Spectro- 
graph (COS) 5 to be installed on the Hub- 
ble Space telescope and the Galaxy Evolution 
Explorer (GALEX) 6 will both be able to de- 
tect Ly— a in the ultraviolet starting in 2003. 
Of course these missions will not be able to 
get spectra for 10 5 quasars, but as we have 
described in §4 we only use the 10 — 12% of 
quasars with logN^j > 19.0 in measuring the 
lensing signal and recovering Thi- Thus if the 
quasars likely to cross this threshold can be 
identified from their Mgll and/or CIV fea- 
tures in a manner identical to that already 
done when looking for DLAS (Rao & Turn- 
shek 2000) it only requires ~ 10 4 spectra to 
measure Nhi in the relevant systems. Since 
GALEX intends to get spectra for 10 4 quasars 
it is even feasible that almost the entire Sloan 
sample can be used in measuring the value of 
Thi- Because the method of targeting Mgll 
absorbers is the most effective way to detect 
DLAS with z < 2, it is likely that a significant 
fraction of telescope time will be spent on such 
an endeavor. About 5 — 10% of all quasars are 
expected to have at least one absorption sys- 
tem that is above our upper bound on Nhi 
(5 x 10 19 cm~ 2 ) and therefore would not be 
included in the lensing analysis. Of Mgll ab- 
sorbers « 50% are above this upper bound 

5 http://cos. colorado.edu/cos 
6 http://www. srl.caltech.edu/galex 
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but the rest will be exactly the higher column 
density LLS that are needed for the lensing 
analysis. Thus the prospects of having UV 
spectra for a large number of quasars with 
19.0 < logN^j < 19.7 is reasonably good. 

One possible impediment to the method 
we have described in this paper would be the 
presence of dust in LLS. Metallicity is usually 
believed to be connected with the dust con- 
tent in galactic systems. If in spite of their 
high ionization states, the LLS do contain 
dust, then dust absorption may be propor- 
tional to Nhi and may act as a counter effect 
to the gravitational magnification. A straight- 
forward way to probe whether absorption by 
dust is important would be to apply the recov- 
ery method from photometry at different band 
filters. Similar values for Thi from analyses in 
different bands should verify the irrelevance 
of dust reddening. If, however, a filter de- 
pendent result is obtained for Thi values, this 
might be an extension of the method which 
would constrain the dust content in LLS. 

6. Conclusions 

We have introduced a new statistical 
method to recover the mass to neutral hy- 
drogen ratio in Lyman limit systems. The 
method makes use of the connection between 
the measured column density of individual ab- 
sorption systems and their contribution to the 
magnification of their background quasars. 

The main effect of the gravitational mag- 
nification is the change of slope at the high 
column density end of the column density dis- 
tribution function. This can be also seen as 
a shift in the quasar luminosity function, a 
shift which depends on column densities. Ac- 
companied effects include a distortion of the 
redshift distribution for the LLS and a small 
shift in the quasars redshift distribution for 
quasars with high column density of total ab- 
sorption. 



The relevant variable for these effects is the 
total column density in LLS along the line of 
sight to a quasar. We have simulated in detail 
the quasar population along with their asso- 
ciated LLS. We have mimicked the observa- 
tional procedure and demonstrated the abil- 
ity to recover the right value of the mass to 
neutral gas ratio, Thi- If Thi is a univer- 
sal constant in these systems then the 95% 
confidence level of the recovered Thi value is 
given by ±140 (y^ n^s'o ) f° r a sample size in 
the range 10 4 < nQso < 10 5 , independent of 
the value of T HI ! For a stochastic Thi with a 
Gaussian distribution with standard deviation 
(7x m , the 95% confidence levels increase only 
slightly because the large number of quasars 
diminishes the effects of stochastic processes. 
A small bias is also introduced because ob- 
served higher values of Thi are preferentially 
brought into the sample and because the mean 
of an unbiased sample is also slightly higher 
since Thi is always positive. However this 
bias is relatively small; for the large scatter 
of o"x HI = Thi the bias is only half as large as 
the 95% confidence level. 

A redshift dependence of the Thi value 
due, e.g., to evolution in the ionization state 
of the LLS can be recovered to a certain ex- 
tent by assuming the functional form of the 
dependence and deriving the free parameters 
of that form. For a reasonable example of a 
redshift dependent model (Haardt & Madau 
1996), the maximum ionization redshift can 
be recovered with an accuracy of ±0.7 These 
scaling relations for the recovery ability are all 
achievable in light of the upcoming quasars 
survey of the SDSS. Follow up observations 
from space should then be undertaken in or- 
der to fill in the low redshift range of LLS, 
unaccessible from ground. 

So far the mass to neutral gas ratio in LLS 
is unknown to more than an order of mag- 
nitude. The method we presented here en- 
ables, for the first time, to constrain this value 
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tightly, hopefully leading to a better under- 
standing of the relation between gas and dark 
matter in the universe. 
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